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Abstract 

Recently, the spectral expansion of finite temperature two-point functions in in- 
tegrable quantum field theories was constructed using a finite volume regularization 
technique and the application of multidimensional residues. In the present work, the 
original calculation is revisited. By clarifying some details in the residue evaluations, 
we find and correct some inaccuracies of the previous result. The final result for con- 
tributions involving no more than two particles in the intermediate states is presented. 
The result is verified by proving a symmetry property which follows from the general 
structure of the spectral expansion, and also by numerical comparison to the discrete 
finite volume spectral sum. A further consistency check is performed by showing that 
the expansion satisfies the cluster property up to the order of the evaluation. 

1 Introduction 



Correlation functions play a central role in the formulation of many-body quantum systems. 
Integrable models presents a unique opportunity to study strongly correlated quantum sys- 
tems in situations where conventional methods break down. Recent experimental advances 
resulted in renewed interest in integrable models, since it is now possible to realize cer- 
tain models with the help of optical and magnetic traps [1, 2, 3, 4] or in low-dimensional 
magnets [5, 6]. 

In a recent paper [7] finite temperature (i.e. thermal) two-point correlation functions 
were constructed using the exact form factors in 1+1 dimensional integrable models. In 
an integrable quantum field theory, the basic object is the factorized S-matrix [8, 9]. The 
matrix elements of the local operators (form factors) satisfy a certain set of equations (the 
form factor bootstrap equations) which follow from general field theoretical arguments 
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supplemented with the special analytic properties of the S-matrix [10, 11, 12, 13]. Solv- 
ing these equations gives the form factor functions, which can then be used to construct 
correlation functions by expanding in the basis formed by the infinite volume asymptotic 
scattering states. 

The form factor expansion of zero-temperature correlations in integrable QFT is very 
well understood. In general, the series has very good convergence properties in massive 
models and can be evaluated numerically to any desired precision [12, 14]. However, the 
problem of thermal correlation functions is much more complicated and has been the sub- 
ject of active research in the last two decades [15, 16, 17, 18, 19, 20, 21, 22, 23]. The form 
factor construction of the spectral series is plagued with problems due to the presence of 
disconnected terms in the expansion, which lead to formally divergent expressions. Fol- 
lowing Balog it can be shown that the divergent parts cancel with contributions from the 
partition function [24]. Nevertheless, it is a very non-trivial task to obtain the correct finite 
answer. Leclair and Mussardo conjectured an answer for the spectral expansion for one- 
point and two-point functions in terms of form factors dressed by appropriate occupation 
number factors containing the pseudo-energy function from the thermodynamical Bethe 
Ansatz [16]. Their proposal for the two-point function was questioned by Saleur [17]; on 
the other hand, in the same paper he also gave a proof of the Leclair-Mussardo formula for 
one-point functions provided the operator considered is the density of some local conserved 
charge. By comparison to an alternative proposal [19], it was also shown that the results 
obtained by naive regularization are ambiguous [20]. 

The idea behind our approach is to use a finite volume setting to regularize the diver- 
gences. In [25, 26] this was applied to one-point functions giving a confirmation of the 
Leclair-Mussardo formula to third order; later a derivation to all orders was also obtained 
[27] . The crucial point is that finite volume is not an ad hoc, but a physical regulator (since 
physically realizable systems are always of finite size), therefore one is virtually guaranteed 
to obtain the correct result when taking the infinite volume limit. The natural small pa- 
rameter for the finite temperature expansion is the Boltzmann-factor e"*"/^ where m is the 
mass gap (which is assumed to be nonzero). The result is an integral series, where the iVth 
term represents iV-particle processes over the Fock-vacuum. The contributions with a low 
number of particles can be interpreted as disconnected terms of matrix elements calculated 
in a thermal state with a large number of particles [23, 27]. In this sense the approach is 
similar to the one used in algebraic Bethe Ansatz [28, 29]. 

Besides correlation functions, the finite volume regularization can also be applied to nu- 
merous other problems. The finite volume form factor approach was extended to boundary 
operators as well [30], which was used to compute finite temperature one-point functions 
of boundary operators [31]. Another application of the bulk finite volume form factors is 
the construction of one-point functions of bulk operators on a finite interval [32] . It was 
also used [33] to construct the form factor perturbation expansion in non-integrable field 
theories (originally proposed by Delfino et al. [34]) beyond the leading order. It also turned 
out that this approach can be applied to quenches in field theory [32, 35, 36]. 

Regarding thermal two-point functions, the finite volume regularization method was 
first applied by Essler and Konik [37, 38]; however, their methods do not have any obvious 
extension to higher order. Despite this shortcoming, their results are very useful as shown 
by their relevance to inelastic neutron scattering experiments [6]. An independent early 
calculation of the one-particle-one-particle contribution can also be found in [39]. 

In [7] we developed a systematic method to compute the finite temperature form factor 
expansion to arbitrary orders. It turns out that the machinery of multidimensional residues 
provides an appropriate formalism to evaluate higher order corrections systematically, and 
this was demonstrated to all orders which involve only intermediate states with at most two 
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particles. To verify the result, we applied two consistency checks. The first of them was that 
the correlator should have a finite limit when the volume is taken to infinity, therefore all 
terms containing positive powers of the volume had to cancel, which was indeed true. The 
second one took into account that for some contributions there are two independent ways 
to arrive at the answer, and agreement between them also provides validating evidence. 
However, for the term which contains the contributions when both intermediate states 
involved in the spectral sum contain two particles, the second one is not available and the 
first is insufficient to check the structure of the result in detail. 

Therefore we decided to provide a numerical evidence, especially since the analytic 
manipulations themselves are rather tedious and complicated, with many possible sources 
of mistakes. As it turned out, the result for D22 reported in [7] is unfortunately incorrect. 
By further investigation, it turned out that some fine details of the residue calculation 
needed to be carried out more carefully. 

Here we report the correct version of the computation and its result, and present the 
final formula for the finite temperature two-point function including all contributions with 
at most two-particle intermediate states. To be confident in our results, we perform several 
checks. First we check that the result for D22 satisfies a particular symmetry property 
following from the general form of the spectral expansion. Then we apply a detailed 
numerical verification of our analytic manipulations, and also verify that the final result 
satisfies the physically required cluster property. 

The layout of the paper is as follows. Section 2 introduces the thermal two-point func- 
tion and the idea of finite volume regularization. In section 3, the methods to evaluate the 
resulting spectral expansion is presented. We re-derive the results of [7] , including the cor- 
rect form of the term D22 and present the full formula of the two-point function including 
all two-particle contributions. The numerical verification of D22 is performed in Section 
4, together with a similar verification for D12 in order to establish a benchmark point for 
numerical accuracy. In Section 5 we prove that the resulting expansion satisfies the cluster 
property, and in Section 6 the conclusions are presented. There are also three appendices: 
Appendix A contains the mathematical formulas used for evaluating the residue contribu- 
tions, while Appendix B contains the end results of the residue evaluations, which are also 
necessary for the numerical comparison. The proof of the symmetry property of the D22 
contribution is given in Appendix C. 

2 Finite volume regularization 
2.1 The thermal two-point function 

A field theory with finite temperature T can be defined using a compact Euclidean (Mat- 
subara) time t: 



A naive spectral sum leads to an ill-defined expression due to the presence of disconnected 
contributions (cf. e.g. the discussion in [7]). However, one can put the system in a finite 
spatial volume L with periodic boundary conditions 



t = t + R where R=l/T 
We are interested in the two-point function in 1 + 1 dimensional field theories: 





(2.2) 



(2.3) 
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so that 



where Tr^, denotes the trace over the finite- volume states, Hi is the Hamiltonian in volume 
L. This expression can be expanded inserting two complete sets of states 

III (e-^^^Oi(x,t)O2(0)) = ^e-^^"(^)(n|Oi(x,t)|m)i(m|O2(0)|n)L (2.5) 

m,n 

where the matrix elements of local operators are also taken in the finite volume system. 
To evaluate it, we need an expression for form factors in finite volume. 

2.2 The form factor bootstrap 

In a 1 + 1 dimensional field theory, the energy and the momentum of an on-shell particle 
is parametrized by the rapidity variable as E = m cosh 9 and p = m sinh 6. For the sake of 
simplicity let us suppose that the spectrum of the model consists of a single particle mass 
m. Incoming and outgoing asymptotic states are defined as: 



\9i, . . . , On) 



|6'i, . . . , OnY"^ : 9i > 62 > ■ ■ ■ > 9n 
|6'i, . . . , 6n)°^^ : 61 < 62 < ■ • • < On 
Integrability leads to factorized scattering, which can be summarized by the relation 
|6'i, . . . ,9k, Ok+i, . . . , 9n) = S{9k — ^A,.+i)|6'i, . . . , 9k+i, 9^, . . . , 9n) 

where S denotes the two-particle amplitude; from this any multi-particle scattering process 
can be obtained by reordering the particles. States are normalized as: 

(^'1^) = 2tx5{9' - 9) (2.6) 

The form factors of a local operator 0{t, x) are defined as 

fK^;, ...,9'J9,,..., 9n)={9[, C|O(0, 0)1^1, . . . , 9n) (2.7) 

With the help of the crossing relations 

FmniO'l, • • • , . . .,9n) = -F^-ln+l(^l: • • • > ^m-ll^m + ^i, . . . , 9n) (2.8) 

n k~l 

+ 2^^(C - Ou) n S{9i - 9k)FZ.,n-A0[, C-il^i, • • • , ^.-1, ^..+1 ■■■.On) 

k=l 1=1 

all form factors can be expressed in terms of the elementary form factors 

F^{9,,...,9^) = {0\OiO,0)\9,,...,9r,) 
which satisfy the form factor bootstrap equations [10, 40, 11] 

Lorentz symmetry: F^{9i + A,92 + A, . . . ,9n + A) = (2.9) 

exp{soA)F^{9i,92,...,9n) 

Exchange: F^{9i, ...,9k, • • • , ^„) = (2.10) 
S{9k — 9k+i)F^{9i, . . . , 9k+i,9k, . . .,9n) 

Cyclic property: (^i + 2t7T, 92,..., 9n) = F^{92, 9^, 9,) (2.11) 

Kinematical poles: —i Res -F^2(^ + 0',0i, . . . , On) = (2-12) 

' n \ 

l-l[SiO'-9k)\F^i9,,...,9n) 

k=l J 

where sq denotes the Lorentz spin of the operator O. There is also a further equation 
related to bound states which we do not need in the sequel. 
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2.3 Form factors in finite volume 

A formalism that gives the exact quantum form factors to all orders in was introduced 
in [25, 26]. The finite volume multi-particle states can be denoted 

|{/l,...,/n})L 

where the are momentum quantum numbers, ordered as /i >■■■>/„ by convention. 
The corresponding energy levels are determined by the Bethe-Yang equations 

l^k 

Defining the two-particle phase shift 5{9) by the relation 

S{e) = -e'^(^) (2.13) 
The derivative of 5 will be denoted by 

due to unitarity, 5 is an odd and ip is an even function. We can write 

Qk{h....A) = mLsm\i~ek + Y^Kh-~ei) = 2Txh , A; = l,...,n (2.15) 

where the quantum numbers take integer/half- integer values for odd/even numbers of 
particles respectively. Eqns. (2.15) must be solved with respect to the particle rapidities 
^fc, where the energy (relative to the finite volume vacuum state) can be computed as 



m cosh 9k (2-16) 



k=l 



up to corrections which decay exponentially with L. The density of n-particle states in 
rapidity space can be calculated as 

p(^i,...,^0 = det:r(") , = , A:,/ = l,...,n (2.17) 

The finite volume behavior of local matrix elements can be given as [25] 

(«,..., 4}|O(0,0)|{/i,...,/„}), = F^^U^'rn + ^^.■■■A + ^^.0..■■■.On) 

+ 0(e-^^) (2.18) 

where 9k {O'k) are the solutions of the Bethe-Yang equations (2.15) corresponding to the 
state with the specified quantum numbers /i, . . . , {![,..., I'^) at the given volume L. 
The above relation is valid provided there are no disconnected terms i.e. the left and the 
right states do not contain particles with the same rapidity, i.e. the sets jfi'i, . . . , and 

^^9'i, . . . , 6'^! are disjoint. 

It is easy to see that in the presence of nontrivial scattering there are only two cases 
when exact equality of (at least some of) the rapidities can occur [26]: 
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1. The two states are identical, i.e. n = m and 

{-^1; • • • 5 I'm} = Ul, ■■■,In} 

in which case the corresponding diagonal matrix element can be written as a sum 
over all bipartite divisions of the set of the n particles involved (including the trivial 
ones when A is the empty set or the complete set {1, . . . , n}) 

({/, . . . i„}mh . . . /„}>. - ^--'"-'f;;'-^";;-"^^^'- + o(e-^) 

p({l,...,n})i 

(2.19) 

where 

is the r-particle Bethe-Yang Jacobi determinant (2.17) involving only the r-element 
subset 1 < ki < ■ ■ ■ < kr < n of the n particles, and 

J^{{k,,...,kr})L = Fl{9k,,...,9kr) 

e-S>0 

is the so-called symmetric evaluation of diagonal multi-particle matrix elements. 

2. Both states are parity symmetric states in the spin zero sector, i.e. 

{/i,...,/4 = {-/„,..., -ji} 

Furthermore, both states must contain one (or possibly more, in a theory with more 
than one species) particle of zero quantum number. Writing m = 2k+l and n = 21 + 1 
and defining 

J^k,i{0[, . . . , d'kl^i, ■ ■ ■ = 

lim F2fc+2;+2(«7r + 6[ + e, . . . ,iiT + 6')^ + e^in - + e, . . . ^in - 6[ + e, 

^n + e,0,9^,...,9l,-9l,...,-9^) (2.20) 

the formula for the finite-volume matrix element takes the form 

...,r„0, -4, . . . , -l[}\0\{h, ...,/,, 0, -Ii, . . . , -h})L (2.21) 
= {P2k+ii9[, . . . , 0, -6^^, . . . , -9[)p2i+i{9i, . . . ,9i,0, -9i, . . . , -9i)j 
X Fk,i{^'i, • • • , 9'i,\9i, . . . ,9i) 

+mL F2k+2i{iT^ + 9'i,...,i'K + 9'^,i'K - 9'^., . . . ,i7i - 9[,9i, ...,9i, -9i, -9i) 
+0(e-^^) 

2.4 The form factor expansion using finite volume regularization 

Using the finite volume description introduced in subsection 2.3 we can write 

(0i(x,t)O2(O))f =i5^C;vAf (2.22) 

N,M 
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where 

CnM = E E {{h-.-lN}\Om\{Jl---JM})LX 
Ii...In Ji---Jm 

{{J^ . . . Jm}|O2(0)|{/i . . . /^})^e^(^-^^)^e-^^(^-*)e-^^* (2.23) 

and Ei^2 and Pi^2 are the total energies and momenta of the multi-particle states . . . In}) l 
and |{ Ji . . . Jm})l- The task is to calculate the sum in finite volume and then take the 
limit L — 7- oo. 

First we classify the contributions into different multi-particle orders following the pro- 
cedure in [38, 7]. Introducing two auxiliary variables u and v (at the end both will be set 
to 1): 

(Oi(x,t)O2(0))f = I E^"^^"'^^^ (2.24) 

N,M 

Similarly for the partition function 

N 

with Zn denoting the A^-particle contribution. The inverse of the partition function is 
expanded as 

N 

where 

Zq = 1 Z\ = —Z\ Z2 = z^ — Z2 
Putting this together we can rewrite the expansion as 

(Oi(x,t)O2(0))f = Y^u^'v^'Dnm (2.25) 

with 

Dnm = '^^Cn-i,m-iZi (2.26) 
I 

The first few nontrivial terms are given by 

Dim = CiM — ZiCq M~i 



D2M — C2M — ZiCi^M-l + (^1 — ^2)Co,Af-2 



(2.27) 



In this way we produce a double series expansions in powers of the variables e~"** and 
g-m(_R-i)_ gince these variables are independent, each quantity Dj^m niust have a well- 
defined L — > 00 limit which we denote as 

Dnm = lim Dnm (2.28) 

and we obtain that 

(Oi(x,t)O2(0))^ = lim (Oi(a;,t)O2(0))f = V Z^^m (2.29) 

N,M 

A similar reordering was also used for the expansion of the one-point function in powers of 
g-mH ^26], and for the boundary one-point function in [31]. It is evident from (2.23) that 
the -Datjv/ with N > M can be obtained from those with N < M after a trivial exchange 
of t with R — t, X with — x and Oi with 02- 
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3 The spectral expansion for finite temperature corre- 
lators 



To evaluate the finite temperature two-point function, it is necessary to evaluate the sum- 
mation over two sets of intermediate states. For a given Cnm this involves an and an 
M particle state. One can start with any of these; to simplify the calculations, it is best 
to start with the one containing the smallest number or particles, and do the other later. 
On the other hand, doing the calculation in the reverse order allows one to cross-check the 
result [7]. 

To evaluate the first summation, a systematic method was given in [7] based on a 
multidimensional residue method. Once this is done, all the singularities from the form 
factors are tamed, and the second summation can be performed by a simple transition 
from the discrete sum to an integral using the density of states. Then, after assembling 
Dnm using the lower Cn'ap coefficients as in (2.26), and taking the limit L — )■ oo the final 
formula for the contribution -Da^m can be obtained. Another quick validity check of the 
calculation is provided by the existence of the infinite volume limit. 

3.1 Converting sums to contour integrals 

For sums over one-particles states |{/})/^ with quantum number / G Z we can substitute 



where 



27re*«iW - 1 



Qi{9)= mL sinh 9 pi (9) = Q[ (9) = mL cosh 9 
and Cj are small closed curves surrounding the solution of 

Qi{9) = 2nl 

in the complex 9 plane. 

For two-particle sums over two-particle states 12})^ with quantum numbers Ji, /2 € 
Z -|- I we can use the multidimensional generalization of the residue theorem to write 



E-E 



d9id92 P2i9i,92 



— — I ir 27r 27r (e'^iiOiM + 1) (QiQ2{di,d2) + 1) 

h>h h>l2 •"-''ih ^ ' ^ ^ 

where C/^/j is a multi-contour (a direct product of two curves in the variables 9i and 6*2) 
surrounding the solution of 

Q^{9i,92) = mL sinh ^1 + 5(^1- ^2) =27r/i 
Q2{9i,92) = mLsmh92 + S{92-9i) =27il2 

where due to the definition (2.13) Ji and I2 take half-integer values, and 
/ dQi dQi \ 

P2(^i,^2) = det ^ ^ \ =m^L'^cosh9icosh92 + mL{cosh9i + cosh92)^{9i-92) 

Since form factors vanish when any two of their arguments coincide, we can extend the 
sum by adding the diagonal: 



h>l2 \h,l2 h=l2/ 
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Figure 3.1: Contour deformation procedure. The black dot shows a singularity not enclosed 
inside the contours following from the spectral sum. 



In the next step, the contours are joined together and opened into straight lines, to a 
product contour whose components in each variable enclose the real axis. However, this can 
only be done by including other poles (apart from the ones needed for the state summations) 
in the interior, which come from singularities of the Q-dependent denominators and of the 
form factors. These must be classified and subtracted. This procedure was discussed in 
some detail in [7], and for one complex variable it is illustrated in fig. 3.1 (for more complex 
variable it must be performed in each variables separately). We shall only outline it for the 
case of the D22 contribution, because of the corrections we make to the previous calculation 
performed in that paper. 

3.2 The D22 contribution revisited 

The D22 contribution is given by 



D22 = lim [C22 - ZiCii + (Z^ - Z2) Coo] = lim 

where 



C22 ^ ZiDii — Z2C[ 



00 



= E E(i^l'^2}|Ol(O)|{Jl,J2})L({Jl,J2}|02(O)|{/l,/2})L<?(^l,^2|^;,^^: 

h>l2 Jl>J2 

with the notation 



K^^^ {{} {} I"!?' "i?' ) g«m3::(sinhi?i+sinhi?2— sinhi9j— sinhi?2)g— i)(coshi?i+coshi92)g— »"i(coshi?'j^+coshi? 

(3.1) 

and where [7] 

Z, = mL f ^cosh^^ie-^^™^^"! 
J 27r 

Z2 = l: /7'^^p2(^i,^2)e-"^^(^°^^''^+™^^'^^)-^ f ^pi(^i)e-2-^-^^''^ 
2 JJ In Ztt 2 J ZTT 

Coo={0,){02) 

Dii = II ^^F^' (^1 + ^2) {§2 + 171. ^1) e^'"^^^-^ 

+ [{0,)Fg^ + {02)Fg^] I 
The rapidities are quantized by 



sinh •O2) ^—m^R—t) cosh i9i g— mt cosh i92 



-mR cosh 1? 



9 



(52(^1,^2) =mLsinh^?2 + 5(^?2-^?i) = 27r/2 (3.2) 



and 



Q'^i'd'^,^^) = mLsmh^[ + 5(^9; -^2) = 27rJi 
Q'2{'d[,^'2) = mLsmh^'2 + 6{^'2 - = 2txJ2 

We perform the Ji, J2-sum first and separate it into a diagonal and an off-diagonal piece: 

= i{Ji,J2} = {luh} term) + ' 

Jl>J2 Jl>J2 

because the finite volume form factor expressions are different for the two types of con- 
tributions. In the second term, the prime indicates that the diagonal contributions are 
excluded. 

3.2.1 The diagonal piece 

This calculation is exactly the same as in [7], so we only highlight the main steps. Starting 
from 

= 5^ (/l,/2|Ol(0)|/l,/2)(/l,/2|O2(0)|/l,/2)i^i?(^l, ^21^1,^2) 



Ii>h 



where 



(/l,/2|O(0)|/i,/2) 



Fj^ (^1, ^2) + (Pl(^l) + Pi (^2)) ^2^ + P2(^l, ^2){0) 



with 



i^2^(^) = F2''{^n,0) 
FZ{^u^2) = limFf (^9i + z7r + £,^?2 + ^7r + £,^92,^9i) 

Writing the sum in terms of contour integrals, after opening the contours and performing 
the large L limit the diagonal contribution becomes 

+ (^^1,^92) e-™^(™^h''i+=°^'^'^2) 

1 rr d^d^[cOsh^,+COsh^2]^ 0) ^-™«(co.h..+cosh..) 

2 Jj 27r 2tt cosh ?9i cosh ^^2 2 v . ; 2 v . ; 
- I ^ [Ff ^ (^vr, 0) (O2) + ^ (^TT, 0) (Oi)] e-™^^-^'^^^ 

+ // Tr^TT^^^ cosh^^i \F^' (zvr, 0) (C2) + if' (^TT, 0) g-^^^-^^^i+^^h'^^) 
JJ ZTT zn 

+ Z2C00 (3.3) 
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3.2.2 The non-diagonal part 

In the non-diagonal part, one can use 



(/l,/2|O(0)|Ji,J2) 



to write 



^nondiag ^ ^ ^ ^ 

h>l2 Jl>J2 



P2(^1,^2)P2W,^'2: 



Xi^g)(^l,^2|^'l,^^) 
C*22 



where 



(^22 - 5Z 



M[M'^ J FfHi92 + ^vr,^9i+27r,^9'i,^9gFf^(i9i + ^7r,i92 + ^vr,^9^,^9'J 



X4?(^l,t92|^'l,^^ 



and the prime denotes the omission of the { Ji, J2} = {-^i, h} term. We can substitute 

Jl>J2 Jl,J2 

since the form factors vanish when any two of their rapidity arguments are identical. 

Now we open the contours to encircle the real axis in and 'd'2. However, that brings 
more singularities inside the contour whose contribution must then be subtracted. These 
can be classified as follows: 

1. Spurious QQ-poles. There are two such terms, which come from including the 
poles with Ji^2 = h,2 or Ji^2 = -^2,1- Their contribution vanishes for L — )■ 00 [7], hence 
the term 'spurious'. However, they must be included in the numerical tests, therefore 
we provide their form in eqn. (B.9). Note that the form factors are not singular in 
this case, although their limits in such points are direction dependent. 

2. QF-poles. In this case one of the integrations has a pole from a form factor, and 
the other one from a Q-term: 

QFI: ^[ = k Q'^i^i,^'^) = 2txJ2 

QFII: d'^ = ^2 ^ Q2(^2,^2) = 27rJ2 

QFIII: = & <5'i(^'i,^i) = 27rJi 

QFIV: = ^2 k g'i(^'i,^2) = 27rJi 

3. FF poles. In this case poles in both integrals come from form factors: 

FFI : ^[ = d'^ = §^ 
FFII: ^\ = 'd'^ = '02 
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The poles of the form factors can be separated by introducing the regular connected part 



where 



- ^1 - ^[ ^2 - ^'2 ^2 - ^'i 
+FZI{^1+ 11^,^2 + i'^W2,'&'i) (3.4) 



A = t{S{A-A)-S{^[-^'2))F^'{A + t7r,^2) 

B = i{S{^[-^'^)S{^2-^i)-l)F^'{^i+i7i,^[) 

C = i{l-S{^2-^i)S{'d[-^'2))F^^{^2 + ni,^'2) 

D = l{S{d[-§'^)-S{d2-^l))F^'{d2+l'K,^\) 

E = i{S{d^-§2)-S{d'^-§\))F^\d2 + i'K,^\) 

F = i{S{^'2-^[)S{^i-^2)-l)F^'{^2 + in,^'^) 

G = Si^i-'d2)Si'&'2-^[))F^^i^i+t7r,&^) 

H = i{S{d'^-§\)- S{d^-§2))F^\d^+i'K,§'^) 

Using the above notation, the pole terms resulting from the form factors can be obtained: 



Ff ^(^^2 + ^vr, + ITT, #1, i?^)Ff ^(t9i + 27r, ^2 + ^vr, = 
F2}A^2 + ^vr, + ni\§\,§'^)FZl{^i + ^vr, ^2 + ni\d'^, d\) 

' E F G 



+Fgl{^2 + ^^,^l+^^^Wl.^'2) 
+FZl{^i + i'K.'d2 + iA^'2.^'i) 



+ 



AH 



(^2 - ^'1) 
AE' 

+ 



^ (^2 
DH 



BG 



'di - §'2 
A 

\_d2 - ^[ 
GF 

' (^i-^'i)' 
AF + GH 



+ 



+ 



B 

^2 - ^'2 ^ 
DE 

(^1 - ^'2? 



^2 - ^'2 
G 



+ 



H 



^2 - ^'i 
D 



AG + BH 



{^2 - ^'1) (A - ^'2 



BE + DG 



+ 



BF + GG 



{^2 - ^'1) (^2 - ^'2 

GE + DF 



{^2 -^'2) {A -^'2) {A -^'2) {A -A) {A -A) {A -^'2) 



(3.5) 



from which one can identify the terms giving QF and FF type singularities. For the residue 
calculation, the formulas of Appendix A can be used. This results in certain differences 
from the result derived in [7], where too simplistic evaluation of residues resulted in some 
inaccuracies in the end result. 

Once the residues are calculated, in the case of the QF terms a further summation 
remains which must be converted into an integral. It has the general form (here written 
for the case QFI): 



E 



G{A,A,A2) 









Q^=27rJ2 
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where d'2 is a solution to 

Q'2(^l,^2) = 27rJ2 

and the case J2 = I2 was omitted since it is a spurious QQ singularity. One can convert 
the J2 summation into integrals using the residue formula 

_ V- /d^G^^9i^^92^ 

2^ f 271 6*^2 (''i .''2) + 1 

Opening the contours and taking care to eliminate the contributions resulting from possible 
poles of the function G lying on the real "d^ axis: 



2ti e*Q2(''i-''2) + 



1 J 271 e'Q'2(^^^^'2) + 1 ^ f 271 e'Q'2(^^''^2) + 1 



^ poles of G /I 

C^2 L^* 



M'2G{^U^2,^'2) G {^1,^2,^2) , j:M'2 Gi^i,^2,^'2 

271 e^Q'2i^^A) + 1 



27X e^'2^^^^^'2) - 


f 1 dQ'^{-di,-d'^) 






a^'2 


1?i,=1?2 



El a(/2 ^ yUl,V2,V2) ,oa\ 
f 271 e'^'2^'^^^^'2) + 1 ^ ^ 

poles o{G p 

"2 

where the second term corrects for the subtraction of the J2 = I2 case and §2 denotes 
the location of the poles of G. The notation ^ corresponds to the straight line contours 
enclosing the real axis as illustrated in fig. 3.1. The full results of the residue calculations 
are given in Appendix B. 

The J2 = I2 term typically is of order 0(1/ L), except for second order pole contribu- 
tions. This results in the following contribution to the QF terms: 



F2^^(^7r,0)F2^^(^7r,0)i^g)(^i,^2|^2,^i) 

[mL cosh 'di — (p — 'f?2)] [ttT'L cosh {^2 — <^ {^2 — '&\] 



[mL cosh '(92 + p {'&2 — "^i)] [mLcosh-i^i + p {-di — '(92)] 



(3.7) 



which is included in QFQ in (B.7). This term was omitted by the calculation performed in 
[7]; its presence is critical for the cluster property. 



3.2.3 Performing the h^h sum and the large volume limit 

We can write 

^nondtag ^ ^22(^1,^2) ^ 1 f _ ) ^22(^1,^2) 



d^9id^92 (722(^1,^2) 

2 ^ y y 27T 27T [e»<3i(''i'^2) + 1] [e«Q2(''i.''2) + 1] 



d^i ^22(^1,^1) Pi(^i) 
2 ^ T 27i [e^Qi(''i'''i) + I]p2(?9i,?9i) 

Since C22 doesn't have any pole we open the contours in the usual way enclosing the real 
axis as illustrated in fig. 3.1. For the L — )■ 00 it is necessary to examine the behavior of 
the Q-functions: 

iQii'&i + ^£1,^2 + i£2) = i^L sinh('i9i + iei) + i6{^i + isi — '^2 — 162) = 

= irriL sinh -^i cos ei — mL cosh 'di sin ei + i6{'di + iei — ^2 — is2) 
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and similarly for Q2 and Q'^ 2- This results in the following limits: 
lim = 

^^00 giQi('9i+jei,i?2+i£2) _|_ X 
1 



lim 

L^oa e«Qi(''i+*ei>'^i+«ei) -|- 1 
1 



lim 7 — 



1, 


ei 


e [0, 


tt] + 


2mT 


0, 




^ r 


zttJ - 


-f In-K 


1, 




e [0, 


tt] + 


- 2mi 


0, 


ei 


e [tt. 


,27l] 


+ 2mi 




Si 


e [0, 


tt] -h 


' 2n7r 


0, 




e [tt. 


,2n] 


+ 2n7r 



Therefore only the upper contours need to be kept, since all other terms vanish exponen- 
tially for large L: 



^^22(^91 + ie, T92 + ie)-l [ ^022(1^1 + ie, A + ie 

2 J 2n 



2 i/ 27r 27r ' ' ' ' 2 J 2ii ' ' ' ' p2{^i + ie.'di + ie) 

In addition, the integrals can be shifted to the real axis. However this leads to singularities 
in the contribution like QF5 (B.6) due to the term containing 

[e^Q'iK,^!) + 1] {^^-^'^) ^ ■ ' 

which can be treated using the identity 

V— =F ^TTo [x) 



X ±i€ X 

3.3 End result for D22 

The terms divergent as L — )■ 00 drop out when including the contribution —ZiDu — Z2CQ0. 
We can also combine some terms by introducing the function 



Ff,(i?i + 27r,^2 + ^vr,^9^,^9;)= (3.9) 

-F2^(^2 + ^vr,/i) 



^{S{^i-^2)-Si^'2-m^o^ 



+ F2,{^l + t7r,^2+^7T\^2,^'l) 
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The end result is 



+ 



X [(1 — S* (^?']^ — di) S {di — §2)) (mxcosh??! — imtsmhdi) 



+ 



^'^F^^ (^2 + ^vr, ^ (^2 + ^vr, ^1) i^ij (^1, ^sl^i, ^1) 

\ i (^1,^2) (02) + Fi^^ (^1,^2) (^1,^21^1,^2) 



+ yy ^^F2^M^vr,0)F2^M^vr,0)J^,^«^(^i,^2|^i,t?2) 

- I ^ (^^,0) (O2) + (^vr,0) (Oi)] 4? 

where V denotes a principal value integral, K^J is defined in (3.1) and Fis is the so-called 
symmetric evaluation of the form factor used in [26] : 

FZ{^i,^2) = limF'^(^9l + ^7^ + e,^92 + ^v^ + e,^92,^9l) 

Note that by introducing F^ss we combined the terms (3.8) into the second integral, hence 
the need for the principal value. 

In (3.10), the underlined pieces are the contributions that are different from the earlier 
calculation performed in [7]. The first underlined term only corrects a typo in [7], where 
this piece was printed with the wrong sign. The second one comes from the subtraction of 
poles of the integrand in (3.6) and the careful evaluation of the principal value term (3.8), 
both of which occur in the manipulation of the QF5 contributions (B.6). 

The third underlined term plays a crucial role in the cluster property. It is the left-over 
from the term 



■m_R(cosh i?i+cosh ^2) 



1 ff d'di d'&o [cosh -|- cosh -do] , n , 

2 ]] 2tx 2tx cosh -i?! cosh ^92 2 v , ; 2 v , ; 

present in the diagonal contribution (3.3), the dependence on 

[cosh T?! + cosh -(92]^ COsh^^i cosh 792 2 

cosh -i?! cosh -(92 cosh^?2 coshi^i 

is simplified by the inclusion of the contribution (3.7), coming from the second order pole 
terms collected in QF6 (B.7). In the large L limit, the terms depending on the cosh ratios 
cancel, leaving us with the last underlined piece in (3.10). As mentioned before, one of the 
mistakes made in the evaluation of D22 in [7] was the omission of this piece. 
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3.4 The full two-point function up to D22 



For completeness, we also give here the lower contributions to the two-point function. 
These are exactly the same as in [39, 7], so we do not give the derivations here. The 
calculations are almost trivial with the exception of D12, where one can use either the 
derivations presented in [7], or follow the steps outlined above, with slight modifications. 
The terms -Da^m with N < M < 2 are 



/ TT'Oi T7' O 2 s'mh-d I— mt cosh I 



27r 

D =1 f '^'^'^ pOi S ) ({^ ^ N g-imx(sinhi?i+sinhi92)-mt(coshi?i+coshi?2) 

02 Q / 2 V 1' 2y \ 2? ly 

Z J ZTT ZTT 



^ imx(sinh i^i — sinh 'd'z) ^—'m,(R—t) cosh i9i ^—mt cosh 1^2 

1 /•/•/• d^i d^l d^'2 ^c.. , , .,Mz.o. 



= 2 ^^^^^ + ^ (vr + e) , ^9;, ^'2) i^s"^^ (^1 + ^{^ + e), #1) 

/AM J„Q/ 
27r 2ti 

(3.11) 

and 1^22 is given in (3.10). In D12 we defined the regular connected form factor function 
Fsrc via the following separation of the kinematical pole terms 

and used the abbreviation 

K^^^ {{} "i?' ) gimx(sinhi9i— sinhi?'j^—sinh 1^2) g—m(i?—t) cosh i?ig—mt(coshi?'j^+cosh 1^2) 

The other contributions Dnm ioi M > N can be obtained from Dmn by exchanging Oi 
with 02 and replacing t ^ R — t, x ^ —x. 

3.5 The symmetry of the D22 term 

The relation between the coefficients -Datm and -Dmat stated above, when applied to D22 
leads to the property that D22 must be symmetric under the following transformation: 

t R-t 

X — !■ —X 
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This is the same as requiring that the result should be independent of which two-particle 
summation is performed first. However, when implementing such a transformation in 
(3.10), the signs of the e terms change, and therefore the contours must be pulled back to 
their original positions. The contour deformation encounters all the singularities on the real 
axis that were treated previously in this section, so the appropriate residue contributions 
must be computed. This computation is relegated to Appendix C, where it is demonstrated 
that the required symmetry property indeed holds, providing the first nontrivial test of the 
result (3.10). 



4 Numerical verification of the analytic results 

The goal of this section is to validate the D22 formula numerically. For this purpose we 
evaluated directly the sum for the two-particle states and compare it with the result of the 
contour integrals. For calculations we used the sinh-Gordon model with the Lagrangian 
density 

1 Tfl^ 

C = -d^<!>d''<i> cosh g'^> (4.1) 

^ 9 

The model contains one massive particle, and its two-particle scattering matrix is simple 
but nontrivial: 

, , sinh 9 — i sin ^ 

S{e) = ^ 

sinh 6 + i sin ^ 

where 

29' 



B 



87r + c/2 

The nontrivial S*- matrix is important, since our formula contains the scattering matrix 
and its derivative in an essential way which we would like to verify. For the fields in the 
correlator, we chose the exponential operators 

e^^s* (4.2) 

normalized to have vacuum expectation value 1, since their form factors are explicitly 
known [41, 42]: 

)(^„ ^„ . . . , e^) = ^ /^'W2...,a:„) j-j ^^^^ _ 



i<j 



i<j 



Xi = e'- (4.3) 

where the polynomials Pn"^ are given by 

^ = [k] 
Pf) = [k]deiM^''\k) n>l 

Mif{k) = [i - j + k]aif_-{xi,X2. . . ,Xn) j = 1, . . . , n - 1 



with 

Hri = ( „, . ' — 1 and \n 



4sin7rP/2\"/' , ^ ^ sin ^ 

—TT—. — and \n\ = b- 

fitir) J ^ ^ sin^ 
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Furthermore, o"!'"'' denotes the elementary symmetric polynomials of n variables defined by 



i=l 



1=1 



ai ' = if / < or / > n 



and the minimal two-particle form factor is given by 



dx 



■ sm 



X 



XitTT — 



2tv 



sinh ^ sinh(l 



f)f sinhf 



sinh X 



(4.4) 



where 



Af = exp 



-4 



dx sinh ^ sinh(l - f ) f sinh | 



X 



sinh X 



(4.5) 



4.1 Evaluating the two-particle sum 

Numerical evaluation of the sum is only possible at finite volume. The factors K^J decrease 
exponentially at large rapidities, so it is possible to choose a rapidity cutoff and restrict 
the summation up to the corresponding Bethe-Yang quantum number. However for large 
volume this quantum number cutoff is still too big and it is practically impossible to 
evaluate the four particle sum. The compromise is to evaluate only the inner two particle 
sum with fixed outer rapidities at moderate volume. This is enough to check the validity 
of all the nontrivial contour deformations and residue manipulation in the £'22 calculation. 
We can write 

QTiandiag _ C*22 (''^l, ''^2) 



with 



^ ^ , ff ■ (tf; + iir, 1 + ix, tfa F?' (tf 1 + «2 + ''T, tf;; 



Jl>J2 



and evaluate C*22 for some given value of -i^i 2, corresponding to a solution of the Bethe-Yang 
equations (3.2) with some quantum numbers {Ji, J2}. 

The parameters for the evaluation can be chosen to help with the convergence of the 
summation, while ensuring that the structure of the expression tested remains general. The 
exponential operators (4.2) in the Sinh-Gordon model can be parametrized by the number 
k that we chose for our evaluations as /ci = 2 for Oi and /c2 = 4 for The essential 
structure of the formula does not depend on this choice. The space-time parameters and 
the temperature were chosen as mx = 0.0, mt = 0.4, and mR = 0.8. Setting mx to zero 
does not hide any important structure of the equation, but makes the expression real and 
that helps in comparing the results with the contour integrals. The sum was evaluated with 
several values for the volume, the sinh-Gordon coupling constant and the outer rapidities: 

ml = (10, 15, 20, 25, 30) 
B = (0.1, 0.2, 0.3, 0.4 ,0.55, 0.7, 0.9) 
rr r. . .5 11 f 11 -51 r-5 -211 (7 -71 fl -1 



2' 2 ' 2 ' 2 ' 2 ' 2 ' 2' 2 ' 2' 2 
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The rapidity cutoff for the quantum numbers included in the sum was chosen as = 
3.0, 4.0, 5.0, 6.0, and the numerical results showed that for the value 'd = 6.0 the discrete 
sum was evaluated within a relative error of less than 10~^^. 

4.2 Evaluating the contour integrals 

To compare the results of the contour integrals with the direct sum, the calculation must 

be performed at the same volumes. In this regime the exponential and power corrections in 

volume L are not negligible, so they must be taken into account. Exponential corrections 

come from the integration on the contours going under the real axis, while power corrections 

come from the total derivative contribution in the second order pole calculation and from 

the { Ji, J2} = {/i, I2} point in the pole and the double integral contributions. The explicit 

formulas can be found in the Appendix B. 

The integration contours run below and above the real axis, and it is important to 

find a choice that is optimal for numerical evaluation. The form factors and hence the 

integrands have poles on the real axis, so it would be better to integrate as far from 

the real axis as possible. However for imaginary parts of rapidities larger than ^ the 
( ft) 

factor Kl";' {^1,^2, ^'1,^2) 

becomes oscillating and exponentially growing in the rapidity 
parameters instead of decaying. Another issue is that the form factors and scattering 
matrices are also evaluated at rapidities that lie out of the physical strip. In the sinh- 
Gordon model the scattering matrix and hence the minimal form factor have poles out of 
the physical strip, with imaginary positions that are proportional to the coupling parameter 
B [41]. Therefore the contour must be chosen to lie between these poles on the one hand 
and the poles on the real axis on the other hand. At the same time it must run as far away 
from all singularities as possible, and also to be closer to the real axis than |. For small B 
this leaves little space for the contours so they run relatively close to the poles, resulting 
in a larger error in the numerical integration. The integration itself was performed using 
Mathematical and the Cuba library for multidimensional numerical integrations [43]. 

4.3 Comparing the results 



B\mL 


10 


15 


20 


25 


30 


0.1 


8.56 X 10-^ 


1.03 X 10-^ 


4.93 X 10-9 


3.29 X 10-^ 


1.31 X 10-^ 


0.2 


7.38 X 10-1° 


3.53 X 10-1° 


6.29 X 10-11 


9.74 X 10-1° 


1.96 X 10-1° 


0.3 


1.74 X 10-1° 


1.13 X 10-1° 


1.25 X 10-1° 


1.26 X 10-1° 


1.26 X 10-1° 


0.4 


1.42 X 10-1° 


1.42 X 10-1° 


1.41 X 10-1° 


1.4 X 10-1° 


1.4 X 10-1° 


0.55 


1.41 X 10-1° 


1.42 X 10-1° 


1.42 X 10-1° 


1.42 X 10-1° 


1.42 X 10-1° 


0.7 


1.33 X 10-1° 


1.34 X 10-1° 


1.33 X 10-1° 


1.33 X 10-1° 


1.33 X 10-1° 


0.9 


1.3 X 10-1° 


1.3 X 10-1° 


1.3 X 10-1° 


1.3 X 10-1° 


1.3 X 10-1° 



Table 4.1: Relative error of the difference between the direct sum and the contour integral 
evaluation of C22 {'di,'d2) with {h.h} = {|, ^} 

Table 4. 1 shows the relative deviation between the direct sum and the contour integral 
evaluation of C*22 (''^i, ''^2) with {/i,J2} = {|)^}- Note that the relative error decreases 
as B grows which can be understood from the conditions for the choice of the integration 
contour mentioned above. Based on the above understanding of the deviations in the 
relative errors for different parameters, and the fact that this pattern of dependence was 

^Wolfram Research Inc, Mathematica, version 8.0, Champaign lUinois, 2010. 
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the same for every value of /i,/2 we checked, it can be inferred that the difference of the 
sum and the contour integration is only due to the numerical errors of integration. 

To provide a further support for this conclusion, the above numerical test was repeated 
for Ci2- The formula of Cu is derived in two independent ways in [7] (depending on 
whether the one-particle or the two-particle summation is performed first), and therefore 
its validity is quite certain even without a numerical test. As in the case of C22, let us 
denote by C12 the result of performing the two-particle summation first with fixed 
rapidity of the one-particle state. Table 4.2 shows the relative deviation between the 
direct sum and the contour integral evaluation of C12 with Ji = 17 as the Bethe-Yang 
quantum number of the one-particle state. The relative deviation has the same pattern 
as for C22 i'&i,^2), and is essentially of the same magnitude. Therefore this evaluation 
gives an independent support for the assertion that the deviations are caused by errors of 
numerical integration. 

As the derivation of D22 from C22 (^1,^2) is almost trivial, the above numerical tests 
also confirm the details of our analjd;ic result for D22- 



B\mL 


10 


15 


20 


25 


30 


0.1 


2.03 X 10-6 


3.56 X 10-6 


3.41 X 10-^ 


1.56 X 10-^ 


8.81 X 10-^ 


0.2 


7.5 X 10-9 


5.03 X 10-9 


1.97 X 10-1° 


1.52 X 10-9 


7.73 X 10-1° 


0.3 


1.48 X 10-1° 


7.76 X 10-11 


3.1 X 10-11 


6.22 X 10-11 


3.13 X 10-11 


0.4 


2.35 X 10-11 


5.29 X 10-11 


5.81 X 10-11 


6.77 X 10-11 


6.59 X 10-11 


0.55 


5.87 X 10-11 


6.61 X 10-11 


6.87 X 10-11 


6.9 X 10-11 


6.96 X 10-11 


0.7 


2.96 X 10-11 


5.32 X 10-11 


6.02 X 10-11 


6.25 X 10-11 


6.38 X 10-11 


0.9 


6.01 X 10-11 


6.48 X 10-11 


6.52 X 10-11 


6.51 X 10-11 


6.52 X 10-11 



Table 4.2: Relative error of the difference between the direct sum and the contour integral 
evaluation of C12 {'di) with Ji = 17 



5 The cluster property of the two-point function 

Another important test of the results is provided by checking that the two-point function 
has the cluster property 

t)O2(0, 0))« ~ (Oi(0, 0))^(O2(0, 0))« 
when the spatial separation x grows large. Using the expansion up to D22 one can write 

(Ol(x,t)(!?2(0,0))« = ^D^A/ 

N,M 

= Doo + Doi + Dw + Du + D02 + D20 + £'12 + D21 + D22 + ... 
For mx ^ 1 the terms containing 

oscillate very fast, and therefore the support of the (multiple) rapidity integrals is restricted 
to the zero measure set 

sinh "dk = sinh 

k I 
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and the integral vanishes. Although this argument looks simple, there is a possible problem. 
Namely, the argument only works if the integrands of the x-dependent terms are all regular. 
A nontrivial example is the term 

in the contribution D12 (cf. eqn. (3.11)), which is in fact regular at d'l = §'2 when the 
two terms inside the braces are added together. In the case of the principal value integral 
in D22 in eqn. (3.10), the regularity of the integrand is ensured by the principal value 
prescription itself. 

As a result, one only needs to examine the terms that are x-independent. We denote 
these by putting a bar over the respective contribution D^m and they read: 



= m{02) 

Doi = ^10 = ^02 = ^20 = 



5n = [(Oi)FfMz7r,0) + (O2)FfM^^,0)]|^e- 



mR cosh 1? 



D12 = ^21 = 

- 1 ff Ml d^2 

J^22 — 



2 // 27r 27r 



(^1,^2) (02)+ (^1,^2) (Oi) 



4s 



^ —mR(cosh i?i+cosh ^2) 



g— 2mi? cosh i9i 



+ // ^^Ff ^ (^vr, 0) ^ (^vr, 0) 
27r 27r 

— [Ff 1 (zTT, 0) {O2) + Ff 2 (zTT, 0) e-2-«™«i^''i 
27r 

The one-point function up to two-particle order is [26]: 

{Of = (0) + |^Ff (^vr,0)e-^-^'^'^^-|^Ff (^vr,0) 

+ - // ^^Ff (^,,^2)e— «(-«i^''i+-^i^''=) + 0(e-3'"«) 
2 // 27r 27r 

As a result one obtains that 

Doo + ^11 + ^22 = m''{02f + O (e"=^™«) 

and therefore the cluster property is satisfied to the given order. Note that the same argu- 
ment shows that the formula for D22 derived in [7] violates the cluster property, providing 
another argument that it needs to be corrected. 



6 Conclusions 

In this work we revisited the finite volume regularization of thermal correlators intro- 
duced in [7], which is based on the finite volume form factor formalism [25, 26]. We have 
shown that the original results for the two-particle-two-particle contribution D22 need to 
be slightly corrected, and presented a modified prescription for evaluation of residue con- 
tributions. As a result, we now have the expansion up to all terms involving intermediate 
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states with no more than two particles. In addition, the result for the nontrivial terms D12 
and D22 was cross-checked with a numerical evaluation of the finite volume summation over 
the intermediate multi-particle states. It was also established that the correlation function 
given by the final formulas (3.10) and (3.11) satisfies the cluster property. In addition, it 
was shown to possess a symmetry property which follows from the general structure of the 
spectral expansion. 

The formalism presented here can be extended to compute any higher correction in 
the series. However, the calculation of D22 is already very tedious, and it is expected to 
become even longer for higher terms. In view of potential applications to condensed matter 
systems, it is likely that the present evaluation would suffice for most of the cases. 

Another reason why the result evaluated up to D22 is interesting is that this is the 
part which generalizes to non-integrable field theories. Breaking integrability in general 
allows inelastic processes; however, below the inelastic threshold the finite volume levels can 
still be described using only the elastic phase shift and the same quantization conditions 
as in (3.2) [44, 45]. Therefore the present computation can be extended self-consistently 
whenever the states dominating the spectral expansion are below the inelastic threshold. 

On the other hand, it remains an open question whether the full expansion for the finite 
temperature two-point function in integrable models can be recast in a form similar to the 
expression conjectured by Leclair and Mussardo [16]. Such an expression would represent 
a partial re-summation of the series, expressing the correlator in terms of Fermi-Dirac 
distribution for the dressed particles corresponding to the representation of the system 
as a free gas of quasi-particles under the thermodynamic Bethe Ansatz [46]. It is known 
that this re-summation is possible for the one-point function [26, 27], even for the case of 
operators located on a system boundary [31]. In that respect, the explicit dependence of 
D22 on the two-particle S-matrix (cf. the first underlined term in eqn. (3.10)) does not 
bode well. The TBA equation contains only the derivative of the phase-shift ip, so any 
dressing term from there is only expected to depend only on ip. Some partial integration 
tricks can be performed to shift this dependence around, but we have found no way of 
eliminating it. On the other hand, it was noticed in [7] that the Di^ contributions allow 
some re-summation by dressing the contributions Do^-i- There is a possibility that a 
redefinition of the form factor terms could help, i.e. if one used some definition for the 
desingularized form factors different from the Frc or Fgg ■ At present it is not known how to 
accomplish this; however, there is still hope for recovering some expression similar to the 
original Leclair-Mussardo conjecture. Evaluation of some higher order corrections could 
shed light on the structure of the series, and the experience gained in the present work 
opens the way to performing these calculations, armed with a numerical method to verify 
the results of the complicated analytic manipulations. We hope to return to this line of 
thought in the near future. 

The present calculation was performed for a theory with a single massive particles. 
Adding more particles to the spectrum is rather straightforward as long as the scattering 
remains diagonal; it is only necessary to add particle species labels in appropriate places. 
For non-diagonal scattering theories, recent progress has made almost all matrix elements 
available in finite volume [47, 48, 49], except for matrix elements involving disconnected 
pieces when the states involved in the matrix element are subject to non-diagonal scatter- 
ing. However, more recently we have solved the issue for two-particle states^ [50], which 
means that the series presented here can be evaluated for general integrable field theories, 
including those with non-diagonal scattering such as the sine-Gordon or 0(3) a models [8]. 

^In fact, the paper [50] presents a conjectured solution for any number of particles, but only the two- 
particle case is backed up by (very strong) numerical evidence. 
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A Residue evaluations 

A.l First order poles 

In this case, the two-dimensional residue formula 

'dzidz2 gizi,Z2) 



g{a,b) 



f J 2m27ri fi{zi,Z2)f2{zi,Z2) det ( ^ 

CaXCt \dZj 

can be applied directly as 

' dzi dz2 



9[Zl,Z2) 



{zi,Z2)={a,b) 



9[a,Z2) 



2-Ki 2m [a - Zi] [eQ^i^^^^^) + i] 



j^dQ2(a,Z2) 




dz2 


Z2=Z* 



CaXCj 

where the (a — Zi) is the pole term coming from the appropriate form factor, and is the 
root of 

^iQ2[a,z*) _|_ ]^ _ g 



A. 2 Second order poles 

These can be evaluated by successive integration, performing first the integral over the 
second order pole coming from the form factor. The fundamental formula to use is 



z=ag{z) g"{a) 3[g"{a)] 



(A.l) 



where for a second order pole g(a) = g'{a) = but g" (a) ^ 0. The terms we need to 
evaluate have the form 



dz\ dz2 



g{zi,Z2) 



CaXCj 

Performing the Zi integral leads to 



2m 2m [a - zi] [e^Q^i^uzi) + i] 



dzo 



dg{zi,Z2) 
dzi 



J 2m I [e*<32{zi,z2) + i] 

Cj 



^ 93(21,22) 




921 


2i=a 


(9Q(zi,zJ) 




9zi 


21 =a 



[eiQ2(2i,22) _|_ X] 

iQ2{a,Z2) ^ dQ(zi,Z2) 
1 



dzo Sy""^ ^2je— z 



zi=a 



2m 



iQ2{zi,Z2) _|_ \Y 



For the second integral, introducing the notation l{zi,Z2) = g{zi, Z2)i^^^^^^ and using 
(A.l) with e^^2(a,22*) ^ _i. 
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dz2 /(a,^2)e*^'^'''^'^ 
27rz [e**32(2i,22) _|_ 1]^ 



d 



dQ2{a,Z2) 
dz2 



dQ2{a,z2) 
dz2 



Putting it together 



CaXCj 



Idzidz2 g{zi,Z2) 

jT 2ni 2ni [a - zif [e'Q^i^^^'^) + 1] 



dg{zi,Z2) 



dQ{zi,z'^) 




dzi 


zi=a 



dzi 



z\=a 



_d_ 

dz2 



'g{a^Z2)i'-^^ 



where again Z2 is the root of 



JQ2{a,zl) 



+ 1 = 



zi=a 



dQ{a,Z2) 

dZ2 



A. 3 Useful formulas for practical evaluation 

Using the form of the Bethe-Yang equations for ^1,^2 

S{^2-^l) = l 



_giQi(i9i,i»2) 

_giQ2()9l,1?2) 



imL sinh '92 



we can easily substitute their solutions into the Bethe-Yang equations for ^'i.,'d'2 



3»Qi(^l,'?2) 

,iQi{^'2,i?^) 



Denoting 



we can also write 



imL sinhi9i 



J,imL sinhi?2 



-5(^1-^2)) 

e-^" ^-^(^2-^2)) 



T (^1,^2,^93 



2-^'l)) = 
5(^91-^2)5(^2-^3 



S{'d2-'dl)S{^l-'d'2) 

Si^2-^i)S{^i-^[) 

S{A-^2)S{^2-^[) 



^-02 A) 
JQ'2i'&[,^2) 



+ 1 



+ 1 



T {^2,^,^2) 
T {^1,^2,^2) 

T(^2,^l,^'l) 

T(^9i,^92,^'i) 



B Pole terms for C22 

Here we list the complete result for the pole term subtractions that appear in 
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(522 (tfl,tf2 



4? (^1,^2K,^'2)X 



-QFl ^^2) - QF2 (^1, ^92) - QF3 ^2) - QF4 (t^i, ^2) 
-gF5 ^?2) - Qi^6 ^92) - FF ^?2) - 2Sgg ^^2) (B.l) 

where the are the contributions from the QF singularities, FF comes from the FF 
singularities and SQQ is the spurious QQ singularity term. As in the main text, the 
notation ^ corresponds to the straight line contours enclosing the real axis as illustrated 
in fig. 3.1. The explicit form of the individual contributions to (B.l) are as follows: 



d^n ( FZli'di + Z7r,^2 + z7r|^ ;,^i)Ff (^2 + it^ ,d'^)K[''J [d^.d^W^^^i] 
27r \^ 



oiQiC^^iA) + 1 



FZli'di + ^vr, ^ 2 + iA^2. A)FZ' (^tt, 0)4? (^1, ^2, ^2, ^1^ 
mL cosh ^^2 + (^2 ~ ^^i) 



(B.2) 



QF2{^,,^2) = (b^ 



( FZl{^^+lT^,^2 + lT^W[.^l)F^\d2 + lTT.^'l)K['t\^l,H^'l.^l, 

27r 



+ 



FZK^I + ZTT, ^ 2 + ^7r|^2, ^l)Ff ^ (ZTT, O)^^!^!, ^2, ^2, ^l) 

niL cosh i?2 + ('*92 — ^i) 



(B.3) 
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X 



— z 



X 



T 271 \ [eiQ'ii^W + 1] - A) 

^21^2,^1) -^2) ->g(^2 -^1)] 

[mL cosh^92 + <^ (?92 - ?9i)] (^^2 - -^i) 
[FfH^9i +27r,i92)Ff2(27r,0) + {Oi ^ O2}] + {^^i ^ ^2}} (B.4) 
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QFm,^,) = I [e-^U^lA) + 1] (^,-^;) 

[mL cosh + {^2 — "^i)] 
X [Ff^(^7r,0)Ff2(i7r,0) + {Oi O O2}] + {^1 ^ ^2}} (B.5) 



njT.M .^d^; / j<^H^i,^2i^;,^i) [^(^; -i^i) - ^(^1-^2)] 



X 



— I 



(^1,^21^2,^1) [S{^2 - ^1) - - ^2)] 
[niL cosh^92 + (?92 - (^^i - ?92) 

[Ff^(27r,0)Ff^(^92 + ^vr,i9i) + {Oi ^ O2}] + {^1 ^ ^2}} 
4?(^i'^2|^i,^9i)5(^i-^92) 

X [Ff (i?2 + z7r,^9i)Ff^(^92 + «vr,?9i) + {Oi ^ O2}] 

+ {^i^^2]} (B.6) 



X 



f ^(^2 + i7r,^;)Ff^(^2 + 27r,^;)j4?(^i,^2|^;,^i) 
27r 1^ [e^Qi('?i-''i) + 1] 

X {T ("(^'i, "(92) [—ma; cosh + imtsvah'di] — mLcosh (-i^i) 

+<^(^9; - di)S{d\ - ^i)s{A - A)} + {A ^ ^2} 

+ ^(ZTT, 0)Ff ^(ZTT, 0)i^i?(^9i, ^921^2, ^1) 

/ [mL cosh -i?! — (yj {-di — ^2)] [mL cosh ~ ^ {^2 — ^1)] 
\ [mL cosh 192 + ^ {'^2 — ^i)] [mL cosh di + (p {di — ^2)] 



+ 



d^; 

I [e^OiW-'^i) + 1] 



[mL cosh + LP - 



X 2 (^^2 + ^vr, i9;)Ff 1 {d2 + ZTT, i?;^ 



+ {^1 ^ ^2} 



+ i^g)(^i,^2|^2,^i)i^2''n^vr,0)Ff^(^7r,0)¥.(^2-^i)' 
1 1 



(B.7) 



[mL cosh^92 + P^ {'&2 — ^1)] [mL cosh^^i + p {i^i — 1)2)]' 
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+ {^91^^92} (B.8) 



SQQ {^1,^2] 



M'^ J Ff 1 (^92 + ^vr, + ztt, i?;, ^9^) Ff ^ + ztt, ^2 + ^vr, ^^'2, ^9;; 



27r 27r 



Xi^'g)(^l,^2|^'l,^'2) ^ (B.9) 



We gave these contributions in their exact finite volume form (i.e. including the full volume 
dependence): albeit they simplify when taking the volume to infinity, and the SQQ term 
does not even contribute in this limit, all terms must be kept in order for the numerical 
verification of Section 4 to work properly. 

C Symmetry of D22 

We want to prove that 1^22 is symmetric under 

t R-t 

X -X (C.l) 

First of all, notice that the diagonal terms contain no t and x factor, and are manifestly 
symmetric under exchanging Oi and O2 ■ So it remains only to treat the non-diagonal 
part. 

C.l The four-integral term 

First we treat the term in (3.10) that contains a fourfold integral. After the transformation 
we change the variables '&i^2 •H- '&'i2^ ^^"^ shift every contour with —2ie. This results in the 
contour now running under the real axis: 

xFf i(?92 + ivr, + in, i}[ - le, d'^ - ie)Ff\di + zvr, §2 + ^vr, d'2 - le, ^[ - ie) 

By shifting the contour above the real axis we can transform this term back to its form in 
(3.10), but during this process we pick up some pole contributions from the poles in eqn. 
(3.5). 

C.1.1 First order pole terms containing F^rc 

Using eqn. (3.5) we can identify a contribution containing Ft^c. In these contribution, all 
poles are of first order, so one can apply the Cauchy formula directly. One set of such 
terms is given by 
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x[5(i?l-^?2)-5(^l-^'l)] 

+ i/ ^^i^c(^2 + ^vr,^9i + z7r|^i,i?'jFf^(^2 + ^7r,^9^)4^^ 

x[5(i9'2-^l)5(^l-^2)-l] 
X [1-^(^1- ^2)5(^2 -^'l)] 

+ 47 ^^4^c(^2 + ^v^,^l + z7^|^2,^2)i^2'''(^l + ^v^,^9^)4? (^1,^921^2,^2) 
X [5(i9'2 - ^2) - ^(^1 - ^2)] 
and a similar contribution from 

XXg)(^l,^2|^'l,^'2) 

= -^^^^Fi^?(^i + .vr,^2 + .vr|^'2,^'i) 

Xirg)(^l,^2|^'l,^'2) 

+ i/ ^^4^c(^i+^vr,^2 + ^vr|^^,^2)Ffn^i+^vr,^^)4^)(^^ 

+ -^j ^Fi^^(^i+m,^2 + ^vr|^2,^;)i^2''n^i+^vr,t?;)4?(^i,^2|^'i,^2) 

x[5(^9;-?92)5(?92-^i)-l] 
+ i/ ^^i^c(^i + ^vr,^2 + ^vr|^^,^i)Ff^(^2 + ^vr,^^)KiJ(^i,^2|^i,^^^ 

X [1-5(^92 -^1)5(^1-^2)] 
+ ^ 1 ^Fi^^(^i + 2vr,^2 + ^vr|^i,^;)Ff^(^2 + ^vr,^;)i^g)(^i,^2|^'i,^i) 

where C± denote integration running above/below from the real axis, in the direction from 
left to right. 

Since the exchange property (2.10) is valid for F^^c in the first and last two variables, 
it can be used to simplify the total contribution to 
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^2-^'l ^2-^2 ^l-^'i ^l-^'2 
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+ 



C 



+ 



D 



^'1-^2 '^2-^2 'd'l-'di ^'2-^1 



28 



1 r J„Q/ [ 



C.1.2 First order poles without F^rc 
The form of these terms is 

Ml d^2 d^'i d^2 



X 



C_ 21T 271 271 271 

AE + DH 



+ 



AF + CH 



AG + BH 



+ 



BE + DG BF + CO 



CE + DF 

Their evaluation is tedious, but straightforward. For example 



AE + DH = - {S{^2 - ^1) - S{^[ - i?'2)) (Si^i - ^2) - S{^'^ - 

X {F^^(^, + in,^'^)F^%^2 + in,^[) + {0' ^O^}} 

and one can evaluate the contributions resulting from the contour shift as 
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d^i d^2 d^?;d?9' 



i^,^,^^(^i,^2i^;,^'2 



AE + DH 

^\d[-d2){^'2-^l) 



271 271 27r 271 

d^i d^2 d^; i^g (^1,^21^1,^1) I 

— — — (SI ' [AE + DH%, ^^^ 



c, 2tc 27X 27\ 



d^i d^2 dd'^Kf^[dxM'&2,'&''. 
27r 27r 27r {d'^ - d^) 



-I \AE + DH%,_,^ 



27r 27r 27r 



d^9i d^92 dd\dd'^.K) 



c. 27r 271 27r 27r 



K,^,^^(^i,^2i^'i,^'2; 



AE + DH 



(^9; - ^92) (^9'2 - ^1) 



d^id^2 /•di9;4?(^i,^2l^'i,^ 



27r 2-K " y 27r (^^'^ - ^2, 
X [5(^92 - ^1) - S{§\ - ^1)] [S{§^ - ^2) - 5(^1 - ^'1)] 



and similarly for the remaining five cases. 
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C.1.3 Second order poles 

We get four contributions which contain a second order pole. All four are the same after 
relabeling the rapidities: 



+ 



Ml M2 d'd[ M'^ 



^. 1.1.1.1 '2.'K 27r 2tt 2tx 



AH 



+ 



BG 



+ 



CF 



+ 



DE 



d^i d§2 M'l ^2 DE 
'c_ 27r 27r 27r 27r {^i - ^2 
ddi (1^2 d^9'i d?9'2 DE 



c+ 2tt 2tt 27r 27r (t^i - ^ 



27r 27r 27r 



^'9=^1 



After performing the differentiation, the final result is: 



Ml M2 M'l M2 DE 



2-K 2-K 2-K 27r (^^i - 



Ml M2 M'l 



F^'{^2 + t7r,^[)F^'{A + t7r,^[)Kl^\^i,A\^[,^i 



X 



2tt 2tx 2ti 

{[mx cosh-i?! — imt sinh-i^i] T {^2-, "di^d'^ T [d'-^^ '&i,'&2) 
[S{^2 - ^i)S{^i - ^'i) - Si^'i - ^i)S{^i - ^2)] <^ {^1 - ^'i)} 



where the notation 
was used. 



C.1.4 Putting together the results 

Putting together the result for the fourfold integral term one obtains 
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1 rrrr d^i d^2 d^'. r.O n a, n/ ^ 



4 JJJJ 2ti 2ti 2tx 2-k 



4 JJJJ 2tx 2ti 2ti 2ti 

^'^^K['^\d,,d2W,,d,)[FZl{§2 + iTi,^i + iA^i^^^^^ 

+ (o^ ^ o^} } - ^i) - s{^2 - M [s{^i - ^2) - s{^i - ^[)] 
.rrd§,dd2 /•d^;4?(^i^^KA)f^o,., ^. ,.^o,(, ^- 



+ {o^^o^]] [s{§2 - ^1) - s{§i - [1 - - §i)s{§i - §2)] 

X {[mx coshi^i — imt sinht^i] T (■(92, i^ij^'i) T '^2) 

+ [S{^2 - ^i)S{^i - - S{^[ - ^,)S{^, - ^2)] {^1 - ^'1)} 

C.2 Other terms 

The following two terms in (3.10) 



2tt 2tt 
ddi f dd[ 



2tx J 2tx 

transform to each other under the symmetry (C.l). 

For the remaining terms in (3.10), after the transformation we redefine ^2 ^ 'd'l to have 
the same K factor as before and get: 
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JJ 2n 2n J 2n 

III U.t/1 U'(/2 -fi 

I 



X 



'jj^^n ^ — ^n^^?^^ — f^^^^ " ^^^^^^^ " ' 

X |Ff ^zvr, 0)Ff ^(^92 + ^tt, #i) + {O, ^ O^} } 



X 



X {r(^?2,^9i,^9;) [-mxcosh^?! +im (/^ - t) sinh^?i] + v?(i92 - ^?i)5(^?2 - - ^^i)} 

C.3 Collecting the terms 

Using the results so far, the transformed D22 is the following: 

+ {F^j-c terms) + (first order pole terms) + (second order pole terms) 



d^^icMs 
27r 27r 



^ (^2 + ^vr, ^i)Ff ^(^2 + ^vr, ^Oi^ff ^ (i^i, t92|^i, ^1) 



+ ^22 
where the F^rc terms are 



+ 



+ /// ^?^^^i?(^i,^2|^;,^i)|F,^i(^i + ^vr,^; + zvr|^2,^i)Ff^(^2 + ^vr,#,) 



27r 27r 27r 



{C^^C^}} (C.2) 
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the first order pole terms are 



— I 



- ' III 1 1 f ^%^{-°' - - . 

+ {0'^0^}] [S{^2 - ^1) - ^(^1 - ^'1)] [1 - - ^i)S{^i - ^2)] 



2tt 271 2tt - 1^1] 



X 



{P^'i'dl + ZTT, ^2)F^'i'd2 + tTT, ^'l) + {Oi ^ O2} 



[S{l92-Vi)S{'&i-l9^}-l\ 



X 



X ■ 



27r 271 ' J 2tt (79; - ?92j 
F^'itn, 0)F^'{A + 27r,i9;) + {Oi ^ O2} } 

{p^'iA + tn, ?9'i)F2^'K + tTT, A) + {Oi ^ O2} } (C.3) 
while the second order pole terms are 

+ Iff ^^^F2^^^2 + -,^;)F2^^(^. + ^v^,/,)<)(^,,^2|^;,^l) 
X {[mx cosh'i?! — imt sinhi?!] T {'d2, 1^1, i^'i) T (i?'^, i^i, '^2) 

+ [S{A - A)S{A - ^'1) - - ^1)^(^1 - ^2)] ^ {A - ^'1)} 



dA dA dA 



F^'iA + i7r,A)F^'{A + i7r,&,)Ki'^\A,A\KA) X (C.4) 



2tt 2tx 2tx 

X {r(^?2,^?i,^?'i) [mxcosh^?! +im {R-t) sinh^Ji] - '^{A - A)S{A - A)S{A - "d'!)} 
C.3.1 Fire terms 

The last two ^4^^ terms (C.2) cancel due to 

FZf{A + iT^.^'i + ^7^|^92,^9l) - i^c'(^2 + n^,A + i^^i.^'i) = 

which can be easily proven from the definition (3.4); the remaining one is combined with 
the first order pole terms below into the F^ss contributions in (3.10). 

C.3. 2 First order pole terms 

The first order pole terms (C.3) can be rearranged into the form 
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+ [O^ ^ O^] } [S{'d\ - ^i)S{^i - ^2) - 1] 

d^,d^2d^; 4?(^i,^2|^i,^;) r^^ +.vr#)F^^f^ +.vr 



using 

S{'d)S{-^) = 1 

These terms are combined with the F/^j.^. terms into the ^4^^ contributions in (3.10). 
C.3.3 Second order pole terms 

The second order pole terms (C.4) have a dependence on i? — t on the last line. One can 
make an integration by parts to transform it into a t-dependence as on the second line. 
Using 

^ ' , ' ^' ' = -mi?sinh(^i)iri? (^1,^21^ ,^1) 
we can make the following transformation in the integrand 

-iT (^92, ^i.^'i) mR sinh ^2\'&'i, ^1) = iT {'^2, ^u^'i) ^"^l^^^^"^^^"^^^ 



(partial integration) -iK[^\di, ^2|^i, ^1) ^ 



and so the contribution is transformed into 

+ III ^^^F,^'i^2 + ^n,^[)F2'^%^. + ^nA)I<5^ 

X {[mxcoshA-^mtsmh^l]T{^[,A.^2) - S{^[ - ^i)S{A - ^2)^ {A - ^'1)} 

With this the second order pole terms (C.4) reduce to the form of the second order pole 
term in (3.10). 

C.4 End result 

After putting together every term, using the definition of F^ss (3.9), and performing some 
simplifications one obtains 

T^trans r-) 

-'^22 — -^22 

which is exactly what was to be proven. 
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